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Abstract 

Motivated by the recent progress of direct search for the productions of stop pair and sbottom 
pair at the LHC, we examine the constraints of the search results on the stop (t1) mass in natural 
SUSY. We first scan the parameter space of natural SUSY in the framework of MSSM, considering 
the constraints from the Higgs mass, B-physics and electroweak precision measurements. Then 
in the allowed parameter space we perform a Monte Carlo simulation for stop pair production 
followed by t; > tX} or i > bx and sbottom pair production followed by 6; > bx? or bi > txi. 
Using the combined results of ATLAS with 20.1 fb™t from the search of £ + jets + Ér, hadronic 
tt + Ér and 2b + Ér, we find that a stop lighter than 600 GeV can be excluded at 95% CL in this 


scenario. 


PACS numbers: 14.65.Ha,14.70.Pw,12.60.Cn 


I. INTRODUCTION 


Observation of a SM-like Higgs boson around 125 GeV by the ATLAS and CMS Collabo- 
rations |1| has had a great impact on the minimal supersymmetric standard model (MSSM), 
in which the naturalness problem is a key issue nin For the electroweak naturalness the 
stops play an important role by cancelling the top quark loop, and thus the LHC search 
for the stop pair production will shed light on this question in the MSSM. The ATLAS 
and CMS collaborations have also made searches for supersymmetric particles with about 
20 fb! of data [94 
SUSY scenario rl 


. The negative results of sparticle search may indicate the natural 


, in which the first two generations of squarks are heavy and the third 
generation of squarks are light. 
In the framework of MSSM, the fine-tuning can be estimated by the Barbieri-Giudice 


measure |3 


A = max{A,, }, 
Olog MZ 
Olog p; | 


(1) 


Pi 


where p; denotes the soft SUSY-breaking terms that enter the minimization conditions of 
the Higgs potential 

M3 _ (mh, + Ea) - (mh, + En) tan? g 

2 p 2) 
with tan 8 = v,/Vvq, Miz, and mir representing the weak scale soft SUSY breaking masses 
of the Higgs fields, u being the higgsino mass parameter, ©, and “gq being the radiative 
corrections to the tree-level Higgs potential. In order to maintain naturalness, the parameters 
LL, My, and X, are expected to be at weak scale, which will impose a constraint on the 


spectrum of MSSM 


m A 
u $ 200 GeV (ST) Vs (3) 


: Mp A 
mi +m? + A? < 600 GeV sin 8 La (=a) E> (4) 
Mhp ) A 


5 kl 
with La = \/3/log(A/TeV) coming from the leading-log approximation of the RGE la, Of 
course, in a given model like the MSSM, CMSSM or natural SUSY, the fine-tuning extent 


2 


depends on the parameters appearing in the left sides of Eqs. (3-5), which can be checked 
by examining the parameter space allowed by current experiments 22 


In this work we focus on natural SUSY and examine the constraints of the LHC search 


results of stop and sbottom_|10, |11]. Note that some studies in this aspect have been 
performed in the literature I|. However, the study in assumed a special scenario 
with A; = 0 and tan = 10 and, as a result, the stop decay is much simplified, while the 
analysis in did not consider the ATLAS sbottom search results. In our study we will 
consider a more general natural SUSY with more experimental constraints. We will first scan 
the parameter space of natural SUSY in the framework of MSSM, considering the constraints 
from the LHC Higgs mass, B-physics and electroweak precision measurements, and then in 
the allowed parameter space we will perform a detailed Monte Carlo simulation for stop 
pair production with tı > t% or tı > by? and sbottom pair production with b= bx? or 
bi > tx; at the LHC, applying the same cuts as in the ATLAS searches. Finally, using the 
combined results of ATLAS with 20.1 fb~! from the search of @+ jets+ Ér, hadronic tt+ Ér 
and 2b + Ér, we will figure out the allowed range for the stop in natural SUSY. 

We organize the paper as follows. In Sec. II, we scan the parameter space and present 
the spectrum of natural SUSY allowed by the Higgs mass, B-physics and electroweak preci- 
sion measurements. In Sec. III we simulate the direct stop/sbottom pair productions and 
examine their constraints on the parameter space of natural SUSY. Finally, we draw our 


conclusion in Sec. IV. 


II. A SCAN OF THE PARAMETER SPACE 


In accordance with the definition of natural SUSY, u should have a value at weak scale 
and the third generation of squarks should not be too far above TeV scale. So we scan these 


parameters in the following ranges 
100 GeV < u < 200 GeV, 100 GeV < (mz, = Mi ,, Mig) < 2 TeV. (6) 


In natural SUSY the gluino mass is supposed to be below several TeV. Since in our study we 
focus on stop/sbottom direct searches, the gluino mass is not a sensitive parameter, which 
is fixed at 2 TeV in our scan. The sleptons and first two generations of squarks in natural 


SUSY are supposed to be heavy, which are all fixed at 5 TeV. For the electroweak gaugino 


masses, we use the grand unification relation Mı : Mj = 1: 2 and take Mı = 1 TeV (in 
natural SUSY M; is usually supposed to be larger than u and the lightest neutralino is 
higgsino-like). Further, we assume A; = A, and Mp = Mg, Which in our scan vary in the 


ranges 
|A| <3 TeV, 1< tang <60, 90 GeV < My, <1 TeV, (7) 


where the lower bound of Ma is indicated from LEP and B-decay constraints. Since the 
component of the stop (t1) will affect its decay modes, we scan the parameter space in two 
scenarios: (I) mz, is lighter than m;,; (II) mz, is lighter than mz,- 


In our scan we consider the following constraints: 


(1) We require the SM-like Higgs mass in the range of 123-127 GeV. We calculate the 
Higgs mass with FeynHiggs2.8.9 and impose the experimental constraints from 
LEP, Tevatron and LHC with HiggsBounds-3.8.0 |26]. 


(2) We require natural SUSY to satisfy various B-physics bounds at 2c level. We use the 
package of Superlso v3.3 |27] to implement the constraints, including B > X,y and 


the latest measurements of B, > uty, Ba > Xt and Bt > rtv. 


(3) We require the ies of the precision electroweak observable such as p, 


sin? Olg, mw and R, [28] to be within the 2c ranges of the experimental values. 


(4) We require the thermal relic density of the lightest neutralino (as the dark matter 
candidate) is below the 20 upper limit of the Planck value 29). We use the code 
MicrOmega v2.4 to calculate the relic density. 


In Fig. [] we show the sparticle spectrum in the parameter space allowed by constraints 
(1-4) in natural SUSY. From the figure we can find the following features and implications. 
Firstly, for both mz, < m;z, and mz, > m;,, the stop (f1) mass is heavier than about 
200 GeV. The lower limit of the stop in case of m;, > m;, is larger than in the case of 
mz, < Mie- The reason is that for mz, > mz, the Higgs mass will receive a negative 
correction proportional to cos 23 log(m;, /m;, ) k and, in order to enhance the Higgs mass 
to 125 GeV, the splitting between m;, and m;, is preferred to be small and a heavier stop 


mass is favored to balance the negative correction. Secondly, the masses of ¥7, ¥? and %8 


are below 200 GeV and are nearly degenerate, which, due to the soft decay products, will 
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FIG. 1: Sparticle mass spectrum in the parameter space allowed by constraints (1-4) in natural 


SUSY. 


weaken the constraints from direct electroweak gauginos searches at the LHC. However, this 
spectrum may strengthen some limits from direct stop/sbottom searches at the LHC since 


the final states b¥y{ (tXī ) have the same signatures as bX > (tX 2). 


III. CONSTRAINTS FROM THE LHC SEARCH OF STOP/SBOTTOM 


Since the decay modes of stop and sbottom will affect the implication of direct searches, 
we first discuss their decay ne ratios. The relevant interactions between stop and 


neutralinos or charginos are given by [32 


Liro = Et fe Pr + SÈ Pa + h.c. ; (8) 
Lage = Gb fE Pr + fRPRXE + he. , (9) 


where P/r = (1 F y5)/2 and 


1 1 : 
i gr Ne liw -+ 370! [Nx] .8 | cos Og — yel Nyo]; iz, Sin Oz, (10) 
2/2 ne x 
E= -=y nle] sin 6; — [Nyo]; 4, cos 67, (11) 
© = plUl cos, (12) 
R = —921V y+] 00867 + Yeas lag sin 6. (13) 
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with y = V2m,/(vsin 3) and y = vV2m,/(v cos 6) being the Yukawa couplings of top 
and bottom quarks, and 6; being the mixing angle between left- and right-handed stops 
(—1/2 < 0z < m/2). The neutralino and chargino mixing matrices N;j, U;; and Vj; are 
defined in The interactions of the sbottom with neutralino or chargino can be obtained 


from above expressions by replacing 0; with 6; and interchanging ys and yp. 
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FIG. 2: The scatter plots of the parameter space allowed by constraints (1-4) in natural SUSY, 


showing the decay branching ratios of stop and sbottom. 


In Fig. [2] we display the numerical results of the decay branching ratios of stop and 
sbottom. From the upper panel we can see that a light stop (t; < 370 GeV) will be 
dominated by decaying to by{. When the stop becomes heavier, the t1 > tXo channels 
are open and the decay fı — bXf is suppressed. For mz, < m;,, Br(t, + bXŤ) is less than 
50%, and Br(t, > t¥?) and Br(t, > tj) are approximately equal which vary between 25% 
and 50%. While for mz, > mgp, Br(t: + bX) is around 50% and Br(t, > tf) can only 


TABLE I: The signals of the ATLAS stop/sbottom direct searches id pi 


Signals of ATLAS stop/sbottom direct searches | Source of each signal in natural SUSY 


pp > it (t > tXo) 
l+ jets + Ér f 


pp > bby (b1 > tX7) 

_ pp > tt, (ty > tx? 

tt(hadronic)+ Ér ea 12) 
pp > bib; (bi > tX7) 


pp— titi (t =>} bxT) 


pp > bibi (bi > 6X99) 


reach 25%. These features can be understood from Eqs. (LOH13) in the chiral limit. From 
those expressions we can see that t; becomes purely left-handed for 6; = 0 and right-handed 


for 6; = +7/2. Since Xi. are higgsino-like, the couplings C, are proportional to y, and 


AX? 2 


nearly equal for both mz, < Mg, and Mmg, > m;,. However, the coupling C};,,+ is dominated 


txt 
by the bottom Yukawa coupling y» for mz, < m;, and by the top Yukawa coupling y: for 
mz, > Mp: So, we can infer that when io tx?» is kinematically allowed, Br(t, + byt) 
will be generally below Br(t; — tX} 2) for mg, < m;, and twice as large as Br(t, > tX 2) 
for mg, > M: For the sbottom decays on the lower panel, they can be understood by 
similar analysis as above. 

In the following, we use the ATLAS direct stop/sbottom search data Lid. kd to constrain 
the parameter space of natural SUSY. The signals of the direct searches and the source 
of each signal in natural SUSY are shown in Table I. In our calculation we simulate the 
signals by MadGraph5 and carry out the parton shower and fast detector simulation 
with PYTHIA and Delphes las We use the anti-k, algorithm to cluster jets and the 
MLM scheme to match our matrix element with parton shower. For the cross section of 
i bb* we use the NLO results calculated by the package Prospino2.1 |38]. Before presenting 
the exclusion limits, in Appendix A we validate our simulations by adopting the same cuts 
flow of ATLAS and comparing our results to the ATLAS results. We can see that our results 
are consistent with ATLAS in the reasonable error ranges. Then, we perform the simulation 
for each survived sample in Fig. 2]and use CL, method to obtain the 95% exclusion limit. 
For each channel i after cuts flow, we can define CL’ as 
_ Pois(n,|b; + s;) 


Li = 14 
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FIG. 3: Same as Fig. [2| but projected on the plane of stop mass versus neutralino mass. The 


samples denoted by bullets (red) are excluded by the ATLAS direct stop/sbottom search data at 


95% CL. 


Here s; is the signal events from our Monte Carlo simulation, and n; and b; are respectively 
the number of the observed and expected SM background events, which can be taken from 
the experimental papers. To account for the systematic uncertainty effect, we allow the 
number of background events in each channel to fluctuate: b; + b;(1 + 6b;). Following the 
analysis in , we convolve the Poisson distribution with a Gaussian weighting factor and 


then integrate over b;: 


where f? is the width of the Gaussian distribution and is treated as the relative systematic 
uncertainty. Besides, we truncate the Gaussian integration in the case of the negative events 


number of background. 


In Fig. 3] we show the constraints of the ATLAS stop/sbottom search data on the parame- 
ter space of natural SUSY. Generally, we can see that the stop mass can be more stringently 
constrained in case of mz, < Mz, than in case of mz, > m;,, which can be understood from 
Fig. 2]and Table 

In the former case, since t + tX{ and b — txy7 are the dominant decay modes for 
moderate stop/sbottom mass (m; > 350 GeV), by using the ATLAS data of hadronic 
tt+ Ér and 0+jets+ Ér, we can almost exclude 350 < mM, < 590 GeV and 350 < mg < 605 
GeV, respectively. But we should mention that since t + ty? and b — tyy have the same 
signal, some samples with a heavier stop can also be excluded by these two channels. But in 
the latter case, the stop mass is hardly bounded by the stop searches due to the suppression 
oft > tXí > and b + tx]. Besides, we note that the sbottom search data of 2b+ Ér can give 
a complementary constraint on the stop mass, which excludes a light stop between 235-350 
GeV in the former case and the exclusion limit can be up to 545 GeV in the latter case. 

In Fig. Alwe combine the above exclusion limits from the three channels with CL, method. 
It can be seen that the data can exclude at 95% CL a stop between 200-670 GeV for 
Mi, < mz, and 250-600 GeV for mg, > m;,. These combined limits are much stronger than 
the limits from single channel. 

Note that since our analysis aim at the natural SUSY where the value of is preferred 
to be small (usually below 200 GeV), our results are relevant for the low u. If u becomes 
large(the lightest neutralino becomes heavy) but the stop mass is fixed, the small mass 
difference between them will make the visible decay products (jets or leptons) of the stop be 
soft. Thus the signal detection efficiency becomes lower so that our direct search limits will 
become weak. We checked and found that when m;, = 500 GeV, u > 300 GeV can escape 
the constraints considered in this work.Besides, we noticed that if u gets heavier, the stop 
should also become heavy to keep the neutralino as LSP, which will suppress the production 
rate of stop pair and escape our constraints. 

The effects of tan @ on our results are shown in Fig[5] in the case of mz, < m;,. From 
this figure we see that the limits get weak for a large tan 3. The reason is that in the case 
of mz, < m;, the decay width of fı + ġib (tı + tX?) is determined by y ~ 1/ cos 
(y ~ 1/sin 8), and thus when tan gets large, the branching ratio of tı —> ty? will be 
suppressed. Then the limits on the stop mass, which dominantly come from the decay 


tı — tx?, will become weak. 


In the case of mz, > Mp, our results are not sensitive to the value of tan 8 because in 
this case the decay widths of tı > {fb and tı — ty? are both determined by y and their 
branching ratios are not sensitive to tan 8. 

In Fig. [6] we display the impact of the direct stop/sbottom searches on the stop sector. 
From the left panel we see that the two regions can be excluded by direct searches at 95% 


CL. The right panel shows that the range of 2 < X;/M,; < 3 can be excluded for mz, < 600 
GeV at 95% CL. 
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FIG. 5: Same as Fig. [4] but showing the dependence on tan £. 
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FIG. 6: Same as Fig. [3| but projected on different planes . 


Since the decay branching ratios of stop will be affected by the parameter A; through the 
mixing angle 6;, we let A; vary in our calculations to examine the its effect on the exclusion 
limit. However, when the stop mass mz, is low (so that the LHC data can set limits), in 
order to get a 125 GeV SM-like Higgs we need a large mz, in addition to a large A+. Then 
from the formula sin(20;) = 2m: Xi/ (mz — mz) where X; = A; — u/ tan 8, we see that the 
value of sin(20;) is small and thus 6; will be close to 0 or +7/2, which means that the stop 
is highly left-handed or right-handed. So the variation of A; can affect the stop mixing, but 
its effect is quite limited due to our requirement of 125 GeV Higgs. We checked numerically 
that the constraints are not sensitive to the value of A; in our case and mainly depend on 


the light stop scale and tan £. 


IV. CONCLUSION 


We studied the constraints of the LHC search results on the stop mass in natural SUSY. 
Considering the constraints from the Higgs mass, B-physics and electroweak precision mea- 
surements, we scanned the parameter space of natural SUSY in the framework of MSSM, 
and then in the allowed parameter space we performed a Monte Carlo simulation for stop 
pair production followed by tı — tX? or tı 4 bYf and sbottom pair production followed by 
bi bX? or bi tx, Using the combined results of ATLAS with 20.1 fb~! from the search 
of L + jets + Ér, hadronic tt + Ér and 2b + Ér, we found that a stop lighter than 600 GeV 


11 


can be excluded at 95% CL. 
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Appendix A: Comparison of our simulation with the ATLAS results 


In the following we present the comparison of our simulation with the ATLAS results. In 
Tables and [IV] we present the results for 2b+ Ér, l+ jets + Ér and hadronic tt + Ér, 


respectively. 


TABLE II: Comparison of ATLAS and our simulation for 2b + Ér when (b1, X?)=(600,1)GeV 


Ér > 150 (GeV) 270 


Pr(j1), Pr(j2), ja (veto) 
Ad¢min(3) 
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TABLE III: Comparison of ATLAS and our simulation for @+ jets+ Ær when (t4, ¥?)=(650,1)GeV 


mss ~- 100 (GeV) 2838 


Ad(jet, E7) > 0.8 2488 
7s > 200 (GeV) 


TABLE IV: Comparison of ATLAS and our simulation for hadronic tt + Ér when 


(41, ¥?)=(600,0)GeV 


e veto 
/# 270.1 303.8 
T > 130(GeV) 


Jet multiplicity and pr 105.1 


pts 30 (GeV) 
Ad(jet, E3) > 1/5 72 82.5 
miss miss,track 
AGED Ep PS) < 0/3 


tau veto and b-tag 33.06 
Sr o o o 


80 (GeV) < m5", < 270 (GeV) 


Er > 150 (GeV) 
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